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1. Introduction 

Random matrix theory (RMT) represents a powerful tool in several statistical problems 
where the relevant degrees of freedom can be encoded as elements of certain ensembles 
of large matrices. Applications cover practically all branches of theoretical physics 
(spectral properties of excited nuclear, atomic and molecular systems, simplicial gravity, 
theory of mesoscopic conductance or Euclidean Quantum Chromodynamics) and are 
more and more frequent in recent interdisciplinary research. Indeed, contemporary 
complex systems are characterized by massive storage of data, usually cast in the 
form of matrices. Typical examples are: huge covariance matrices for financial 
and/or economic indicators, incidence matrices for information/retrieval search engines, 
adjacency matrices for networks, stochastic matrices describing frequencies of alleles in 
various gene populations, frequency response matrices in wireless telecommunication 
systems, to mention only few. Recently, concepts of free random calculus [H [2] allowed 
to build in parallel to Classical Probability Calculus (hereafter CPC) a new framework 
- matrix-valued probability theory, hereafter named Free Probability Calculus (FPC), 
and allowed to recast several mathematical features of large size random matrices in an 
intuitive and user-friendly way. 

Most of the above-mentioned applications of random matrices of free random 
variables correspond to the "static case", when we are only interested in spectral 
properties of certain ensembles, and not in the dynamical properties of evolution of 
such ensembles as a function of some external parameter, representing either real 
time or the length of a disordered wire, temperature or the numerical value of 
the coupling constant. In all these realms matrix ensembles are becoming evolving 
dynamical systems. This is the problem we are addressing in this note: using the 
parallel between standard probability calculus for one-dimensional random variables and 
noncommutative probability calculus of large matrices (FPC) we consider the simplest 
matrix-valued stochastic evolutions - Brownian walks. Note, that in our case the 
stochastic increment is represented by a whole matrix, so we deal rather with Brownian 
crowds than with single walker. 

The formulation of this problem was addressed in our previous papers [3l HJ E], 
where we have also given references to related literature on this subject. Here, 
after very brief familiarization with such concepts of free variables as multiplicative 
S transformation (section 2), we review two simple cases of stochastic multiplicative 
evolutions considered recently in the literature [3l [5] (unitary multiplicative evolution 
and hermitian multiplicative evolution) and we establish a new link between their 
spectral properties. 
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2. S transformation 

One dimensional diffusion process for geometric random walk (in the absence of drift) 
is given by the stochastic differential equation 

— = adx T (1) 

y 

where y is random and dx T represents Wiener measure. By averaging y(r) we recover 
well known solution for the probability density, so-called log-normal law 



p(y,T\yo,0) = , exp 



Qog(y/yo) + Wt 



(2) 



The corresponding analogue in the space of matrices may be written as the product 
of random matrix- valued variables 

K 



i 

where Xi represent large (infinite) hermitian Gaussian free matrices belonging to GUE 
with the unit variance and infinitesimal time step is equal to r/K. 

To solve this multiplicative diffusion process we notice, that since matrices are 
non-commuting objects, a standard exponentiation trick relating algebraic to geometric 
Brownian walk fails. This is why we have to use the concepts of free random variables. 
We proceed in the following steps. 

First, we consider an ensemble built of large N x N matrices M, for which the 
pertinent resolvent (Green's function) is generically given by 

G W =/jM e -"l»)il(-i i ), (4) 

In the case of Gaussian ensemble, potential V(M) equals |TrM 2 . For large z, the 
Green's function is the generating function for spectral moments (we do not consider 
here so-called Levy ensembles, for which moments do not exist). In the Gaussian case 
G(z) = (z — \J z 2 — 4)/2, so the imaginary part yields Wigner's semicircle spectral 
distribution and expansion in 1/z generates all the spectral moments of this distribution. 

Second, note that the product of hermitian matrices is, in general, not a hermitian 
matrix anymore, which means that the eigenvalues start to diffuse over the whole 
complex plane, and do not stay on the real axis. This problem was recently solved by us 
using diagrammatic tools of random matrix theory [3]. Note however, that in the case 
when we are only interested in moments of the resulting non-hermitian distribution, we 
may still use the concepts of the generating functions, i.e. resolvents, and considerably 
reduce the technical complexity of the whole problem. This is due to the fact that even 
for non-hermitian ensembles, whose eigenvalues lie in bounded domains of the complex 
plane, the Green's function (j3J) is analytic around z = oo and thus serves as a generating 
function for the moments. 

In order to carry out this program in practice, we use here a remarkable tool of FPC 
known as multiplicative S transformation, allowing us to find a generating function for all 
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moments of the product, provided the spectral properties of the individual ensembles are 
known. The operational procedure is quite simple: Let us assume that we know Green's 
Gi(z) functions for each of the hermitian ensembles Mj. We construct an auxiliary 
analytical function x% f° r each of them, defined by the condition 

n -l = z (5) 




and then we redefine the result, getting Si(z) = ^-^Xi( z )- The resulting S-transform is 
multiplicative, i.e. 

S lx2 (z) = S 1 (z) ■ S 2 (z) (6) 
and represents the multiplication law for the spectral cumulants. 



Applying the concept to individual matrix-valued entry 1 + yr/KX, we obviously 
recover Wigner's semicircle shifted by 1 and with rescaled variance (3 = r/K. A short 
calculation [I] yields the desired S-transformation for the product 

K 



S Y (z) = lim S K = hm J] Si = hm 1 — (7) 



2tz/K 



The surprisingly simple answer reads 

Sy(z) = exp(— tz) (8) 

Note that despite the fact that we have considered initially only Gaussian ensembles 
(with unit variance), the above result is valid for any hermitian random matrix ensemble 
provided its second moment is finite and normalized to 1. In such a way it represents a 
new central limit theorem for the product of infinitely many hermitian entries. 

Let us now consider a quite different multiplicative evolution, solved recently by 
two of the present authors [5], namely the spectrum of a unitary operator, represented 
by an infinite product of unitary, free (independent) increments. Now the evolution 
operator is represented as 



Y T = Ue^UJ-H i j (9) 

i 

where one assumes only that the second spectral moments of hermitian Hi are finite. 
Note that in this case the spectrum always stays on the unit circle, since the product 
of unitary matrices is unitary. In this case, the S-transformation is again surprisingly 
simple and the result reads [5] 

S(z) = exp(r(^ + 1/2)) (10) 

where we have set the arbitrary variance of the initial distribution to one, to simplify 
the comparison between the unitary and hermitian cases considered here. 



3. Dynamical evolution of both ensembles 



To distinguish between both cases of unitary and hermitian S transformations, we 
denote them by U and H subscripts, respectively. Since Sh{z) = exp(— tz) and 
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Su(z) = expr((z + 1/2), we notice that 

S H (z)Su(z) = e T ' 2 (11) 

i.e. spectral properties of moments are very closely related. To see the relation to 
moments in more detail, we recall that the Green's function G(z) for the unitary diffusive 
matrix ensembles [5] generates moments 

11 00 

G(-)-l = J2™ n z n (12) 



Z Z n=l 



where 

m n (r) =< TrY T n >= ±- I p(0, r)e m6 d6 ( 13) 

271 J 

Since the multiplicative Voiculescu S transformation is related to the functional inverse 
of the generating function: 

J G[-ir]-l = * (14) 



where S(z) = -^x(z), one can apply standard residua trick for inverting the series [6], 
getting that 

(15) 

In the case of Su(z) the integral representation for moments reads: 

mn = -f + l/z) n eM-nr{z + 1/2)) = -e~^ L^rn) (16) 

n J 2m n 

where we used the integral representation of the Laguerre polynomials. This result 

resembles the structure observed in the Gross- Witten type models [7J. A closer look 

shows, that the similarity to Gross- Witten models or string models is actually exact, 

provided that the corresponding coupling constant in this model (or the surface of the 

string in string models) is identified with the evolution parameter r (in units where the 

second spectral moment is one). As far as we now, this simple observation relating the 

generic models to diffusive properties of large matrices was never explicitly stated in the 

literature. 



4. Phase transitions 



It is important to notice that both considered models have phase transitions. In 
the case of unitary diffusions, the phase transition is visible e.g. from the large n 
behavior of moments. In this case, asymptotics of Laguerre polynomials known as 
Plancherel-Rotach limit confirms the critical value r c = 4, above which moments decay 
exponentially and for times smaller than critical, moments oscillate with amplitude 
^-3/2 _ These results were explicitly derived in the context of two-dimensional QCD 
by Gross and Matytsin [8], provided we identify string surface A with r. Resulting 
phase transition is third order (Douglas-Kazakov type), and generic for all types of 
Gross- Witten type models. It corresponds to a geometrical phase transition, when 
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diffusing spectra, originally concentrated at single point or the finite arc of the unit 
circle, spread in both orientations of the circle to meet and to close the gap previously 
free of eigenvalues. Our duality relation ( TTTT) suggests the presence of phase transition 
at r c = 4 also for the case of the the product of hermitian ensembles. This was the case 
observed analytically and confirmed numerically in [3]. The phase transition there was 
of "topological nature" again - it corresponded to the appearance of the hole on the 
plane of diffusing complex eigenvalues. 



5. Burgers equations 



Last but not least, we address the issue of dynamical equations describing above- 
mentioned phenomena. To understand the structure of the a-priori unknown equations, 
we consider first the simple additive evolution of Wigner semicircle. This is the case 
when the length of support averaged spectral distribution of eigenvalues has to grow as 
y/r, in analogy to Einstein-Smoluchowski formula for "classical diffusion". This process 
was already considered in 1962 by Dyson [9], who was considering stochastic differential 
Langevin equations for temporal behavior of the eigenvalues 

dHT) ^ d ^ + ^ jl - ( - j L_ )dt (17) 

where the first term represents Wiener process and the second, originating from the 
measure (van der Monde determinant) represents deterministic, strong repulsive force 
between any pair of eigenvalues. For large N individual Langevin-like equations may 
be represented by the equation for average spectral distributions p(A, r), i.e. analogs of 
Fick/Fokker-Planck equations known in classical case. They read [10] 



d T p(X, t) + d x 







therefore are non-linear. An alternative way to recast these equations is to use the 
Green's function. Since imaginary part of the Green's functions is related to spectral 
density, matrix- valued "diffusion equations" take the form 

d T G(z,T)+G(z,T)d z G(z,r) = (19) 

known also as complex Burgers equation. 

We can return now to the case of the hermitian multiplicative diffusion. Knowing 
Sh{z) = exp(— tz) we read out the Green's function for the ensemble Y 

l + f(z,r) 



G y (z,t) 

z 

where / is the solution of the transcendental equation 

Z f = (i + fyf 



(20) 
(21) 
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Note, that by definition, Gy{z,t) is a holomorphic function (defined for large z) 
generating spectral moments of Y T , i.e. 

OO 1 

G Y (z,r) = j:trY?— T (22) 

n=0 ^ 

Let us remind [4J that Sy fulfills a simple equation 

zd z S - rd T S = (23) 
By differentiating (T2"Tj) with respect to r and z we get the following evolution equation 

d T f + zfdj = (24) 
or equivalent ly 

drf + fdv nz f = Q (25) 

supplemented by the boundary condition /(r = 0, z) = (z — l) -1 . This is the Burgers 
equation where spatial evolution is governed by log z. Note however, that this time 
the semicircle solution of the additive Burgers equation (fl9l) . with z replaced by lnz, 
does not fulfill f l24|) . It should not be puzzling for the reader that similar Burgers-like 
equations describing the evolution of the spectra one encounters also in the case of 
unitary evolutions [5j[8]. 



6. Summary 

In this paper we discussed some problems of Brownian-like walks in the space of the 
spectra of large matrices with the help of powerful tools of free random variables 
calculus. New results presented in this note include: (i) a rather unexpected link 
between the moments of unitary evolutions and hermitian evolutions; (ii) an exact 
correspondence between diffusive character of evolving large unitary matrices and the 
whole class of generic models of Gross- Witten type known in two-dimensional Quantum 
Chromodynamics. We hope that the language of matrix-valued evolution presented 
here does not only form an exciting and rich mathematical structure, but may also help 
to understand basic concepts of "new statistical physics" of large and non-commuting 
objects. 
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